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ABSTRACT 

We have assessed the significance of Tuhy & Verheijen's (1997) bimodal Ursa Ma- 
jor Cluster spiral galaxy near-infrared surface brightness distribution, focussing on 
whether this bimodality is simply an artifact of small number statistics. A Kolmogorov- 
Smirnov style of significance test shows that the total distribution is fairly represented 
by a single-peaked distribution, but that their isolated galaxy subsample (with no sig- 
nificant neighbours within a projected distance of ~ 80 kpc) is bimodal at the 96 per 
cent level. We have also investigated the assumptions underlying the isolated galaxy 
surface brightness distribution, finding that the (often large) inclination corrections 
used in the construction of this distribution reduce the significance of the bimodal- 
ity. We conclude that the Ursa Major Cluster dataset is insufficient to establish the 
presence of a bimodal near-infrared surface brightness distribution: an independent 
sample of ~ 100 isolated, low inclination galaxies is required to establish bimodality 
at the 99 per cent level. 

Key words: galaxies: spiral ~ galaxies: general - galaxies: statistics - galaxies: clus- 
ters: individual (Ursa Major) - galaxies : photometry - galaxies: evolution 



1 INTRODUCTION 

The distribution of spiral galaxy disc surface brightnesses is 
of great observational and theoretical importance in extra- 
galactic astronomy. Observationally, knowledge about the 
present day disc galaxy central surface brightness distri- 
bution (SBD) allows us to better understand the impor- 
tance and implications of surface brightness selection effects, 
which is important for e.g. understanding the local galaxy 
luminosity function (e.g. Dalcanton 1998; de Jong & Lacey 
1998), or in comparing galaxy populations at low and high 
redshifts (e.g. Simard et al. 1999). Theories of galaxy forma- 
tion and evolution link the SBD with the angular momen- 
tum, mass and star formation history of galaxies: indeed, the 
SBD has been used as a constraint on some galaxy evolution 
models (Dalcanton, Spergel & Summers 1997; Mo, Mao & 
White 1998; Baugh et al. 1998). 

Previous studies have generally found a SBD that can 
be described as flat or slowly declining with surface bright- 
ness, modulated by a sharp cut-off at the high surface bright- 
ness end, normally associated with the Freeman (1970) value 
(e.g. de Jong 1996c; McGaugh 1996). 

TuUy and Verheijen (1997; TV hereafter) derived the 
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SBD of a complete sample of 60 spiral discs in the Ursa 
Major Cluster from an imaging study in the optical and 
near-infrared (TuUy et al. 1996). Their data suggested an 
inclination-corrected near-infrared bimodal SBD (as opposed 
to the conventional unimodal distribution described above) 
with a peak associated with the Freeman value, and a sec- 
ond peak ~ 3 mag arcsec"^ fainter, separated by a gap (Fig. 
1, solid line in panel a). This would have far-reaching im- 
plications for our understanding of galaxy formation and 
evolution. A bimodal SBD would imply that there are two 
discrete galaxy populations: a high surface density popula- 
tion and one with surface densities a factor of ~ 10 lower 
(assuming comparable mass to light ratios for the two pop- 
ulations). Correspondingly, the masses, spin distributions 
and star formation histories of these two families of galax- 
ies would adopt one of two preferred states, and plausible 
mechanisms for generating these discrete differences in e.g. 
star formation history, or spin distribution, would have to 
be found. 

Why this bimodal distribution was not found before is 
an important question. TV address this problem at length: 
here we will summarise their main arguments. First, the 
bimodal SBD is based on near-infrared K' data and, in 
fact, is only really apparent there. Near-infrared photom- 
etry is much less susceptible to the effects of dust extinction 
than optical photometry. Furthermore, high surface bright- 
ness galaxies are typically redder than low surface bright- 
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ness galaxies (de Jong 1996b), which accentuates the gap 
in the near-infrared. TV demonstrate that in the optical 
B band the bimodality is washed out by these effects: the 
bimodality in the optical is consistent with the bimodality 
in the near-infrared, but is less significant (accordingly, wc 
do not consider TV's optical SBDs here). A bimodal SBD 
would therefore remain undetected in optical studies (e.g. 
McGaugh 1996, de Jong & Lacey 1998). 

Second, the TV SBD has been determined from a com- 
plete sample of galaxies at a common distance, making se- 
lection effects much easier to understand. This may explain 
why de Jong (1996c), using K band photometry for a sample 
of field spiral galaxies, did not report any sign of bimodality 
in his SBD. 

In an attempt to investigate whether the bimodal dis- 
tribution could be caused by environmental differences, TV 
defined an isolated subsample, consisting of 36 galaxies that 
have a projected distance > 80 kpc to their nearest neigh- 
bour (assuming a distance to the Ursa Major Cluster of 15.5 
Mpc, corresponding to Ho = 85 kms^^ Mpc~^). They found 
that the gap is more pronounced in this isolated subsample, 
and attributed intermediate surface brightnesses found in 
their total sample to interactions (Fig. 1, solid line in panel 
b). 

Unfortunately, TV did not attempt to robustly estimate 
the significance of the apparent bimodality. Their only esti- 
mate is based on the isolated subsample where they compare 
the number of galaxies in the 1.5 mag wide gap (between ~ 
17.5 and ~ 19 magarcsec"^; 3 galaxies) with the number of 
galaxies that they expected to see in that gap (20 galaxies). 
This expectation value was derived from a model unimodal 
distribution normalised to reproduce the peak number of 
galaxies in the high surface brightness bins, not the total 
number of galcixies. Properly normalising the model distri- 
bution lowers the expectation from ~ 20 galaxies to ~ 12 
galaxies. 

In this paper, wc investigate the significance of the 
bimodal near- infrared SBD, focussing on whether this bi- 
modality is simply an artifact of small number statistics. In 
sections 2 and 3, we use a Kolmogorov-Smirnov (KS) style of 
test to address this question, investigating the significance of 
both the total and isolated samples. In section 4, we discuss 
the results, focussing on the role of inclination corrections. 
We present our main conclusions in section 5. 



2 CONSTRUCTING A NULL HYPOTHESIS 

In order to test the significance of the bimodality in the 
SBD, it is necessary to choose a null distribution to test the 
observations against. Ordinarily, it would be best to choose 
some kind of 'best bet' observationally determined or the- 
oretically motivated null hypothesis: failure of this null hy- 
pothesis would indicate that it does not accurately describe 
the dataset. However, in this case, we wish to quantify the 
significance of a particular feature in the observational dis- 
tribution, namely the bimodality in the Ursa Major Cluster 
SBD. For this reason, it is necessary to be careful about 
choosing a null distribution: if a null distribution was cho- 
sen which did not adequately describe the general features 
of the SBD, then a high significance (i.e. a low probability 
of obtaining the data from the null) would not represent a 
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Figure 1 . The SBD for the total sample of 60 galaxies (solid line 
in panel a) and for the isolated subsample of 36 galaxies (solid 
line in panel b). Also shown in panels a and b is a normalised 
representative null distribution (dashed line). Panel c contains 

our rc-fittcd isolated galaxy subsample, and panel d shows the 
low inclination isolated galaxies only. Panel c contains a repre- 
sentative model distribution taking selection effects into account 
(solid line), and shows the same distribution without taking into 
account the selection effects (dashed line, with an arbitrary nor- 
malisation). Panel f shows a 36 galaxy realisation of the single- 
peaked model distribution, chosen to appear reasonably similar 
to the isolated galaxy subsample. Panels g and h show the cumu- 
lative distribution function for the total and isolated subsamples 
respectively. Solid lines depict the observed cumulative SBDs and 
dashed lines represent the model distribution. 



large degree of bimodality, but rather would indicate that 
the null hypothesis was a poor match to the general shape 
of the SBD. 

For this reason, we have chosen to fit a relatively simple 
model distribution to the data (taking into account selection 
effects), allowing us to determine the significance of the bi- 
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modality by minimising the mismatches between the general 
shapes of the two distributions. 

For the unimodal model, we use a slightly modified ver- 
sion of the optical SBD presented in McGaugh (1996): 

logio$(Mif',o) = m(/iK',o - Mk',o) 

where ( rn = ruf ii fiK' ,o > fJ-h ,o 
I m = 2.6 if Hk',0 < Hk',o 

The critical surface brightness fi*j^/ g and the faint-end slope 

m/ are adjustable parameters which are fit below. ^ For 
reasonable values of m/, the number of galaxies per unit 
magarcsec"^ only slowly changes with surface brightness, 
and falls quite sharply to zero above a given 'critical' surface 
brightness fj,^, q (see e.g. the dashed line in panel e of Fig. 
1)- 

To take into account the selection effects that affect 

TV's observations, it is also necessary to adopt a scale 
size distribution: we assume that the scale size distribu- 
tion is constant (per logarithmic interval) over the interval 
logjo = 0.8 to logj^o Hk' = 1.7 (~6.3 arcsec to 50.1 arc- 
sec, or, at an adopted distance of 15.5 Mpc, ~0.5 to 3.8 
kpc). 

Selection criteria are modelled by choosing only galaxies 
with isophotal magnitudes rriK' < 12, and isophotal diame- 
ters Dj(t > 60", both measured at the 22.5 K' mag arcsec"^ 
isophote. Though these criteria differ slightly from those 
used by TV (as their selection was made in B band rather 
than K' band), they are comparable to better than 0.5 mag, 
assuming B — ' ~ 3 for faint late-type galaxies (de Jong 
1996b), and accurately follow the surface brightness and size 
limits on TV's fJ,K',o — logj^g hni plane. We emphasize that 
we do not hope to somehow represent the 'universal' K' 
band SBD with this function, we simply aim to construct a 
plausible, simple SBD that provides a sensible single-peaked 
fit to the observed bimodal Ursa Major SBD. 

In panel e in Fig. 1, we show a 10^ galaxy Monte Carlo 
realisation of a model distribution with q = 16.3 and 
ruf = 0.0 (also in panels a and b, and in cumulative form 
in panels g and h as a dashed line). We also show the same 
model distribution without the selection effects (with arbi- 
trary normalisation) in panel e as a dashed line. 



3 RESULTS 

In the above section, we described the single-peaked distri- 
bution that we will fit to the observed SBD. As we want 
to take into account selection effects in the null distribu- 
tion, we cannot make a straightforward fit, but have to use 
Monte Carlo realisations of the null distribution and asso- 
ciated selection effects. We generate a grid of Monte Carlo 
unimodal distributions, each containing 10'' galaxies (later 
renormalised to contain the same number of galaxies in the 
observed SBD), and each with different values of q and 
mj (with grid stops of 0.05 in /i^/ g and 0.02 in m/). The 
best fit unimodal distribution is determined by minimis- 
ing the Cash statistic C = — 2 ^"^Y" log Pi(n, A), where 

•t McGaugh's original values were /ig = 21.5 Bj magaxcsec"^ 
and rrif = —0.3, which are roughly consistent with the best fit 
results described in section 3. 



Pi(n,\) is the Poissonian probability of observing n galax- 
ies in a bin which the unimodal model predicts should have 
a mean galaxy number of A (Cash 1979). The Cash statistic 
has the advantage that the underlying distribution does not 
have to be Gaussian. As the model grid used to determine 
the best fit is itself a Monte Carlo realisation of the true 
underlying model grid, the best fit parameters will depend 
slightly on the model grid realisation used. This ultimately 
affects the significances that we will derive later. We have 
therefore run the code 10 times, and thus constructed 10 
reahsations of the entire model grid, and derived a best fit 
for each of these realisations. The best fit unimodal distri- 
bution typically has parameters around /i^/ g ~ 16.35 ±0.05 
mag arcsec"^ and a faint end slope m/ ~ 0.05 ± 0.08, where 
these uncertainties are the RMS variation in the fit param- 
eters. 



To estimate the significance of the departures of the 
observed SBD from the best fit model distribution we use the 
KS statistic dmax, defined as the maximum distance between 
the cumulative distributions of the model and the data (i.e. 
the maximum distance between the solid and dashed lines 
in panels g and h of Fig. 1). The dmax.obs between this best 
fit model and observed SBD is then measured. 



However, because we have fitted a single peaked distri- 
bution to the bimodal SBD, it is not possible to convert 
the measured dmax.obs into a probability in the standard 
way (e.g. Press et al. 1986). ^ Therefore, it is necessary to 
use Monte Carlo simulations to convert the measured maxi- 
mum distance into a probability. Accordingly, we generated 
100 random simulated datasets from the best fit model, and 
subjected them to the same fitting and maximum distance 
measuring procedure as the data. In this way, we determined 
what the distribution of maximum distances should be from 
the simulations. The fraction of simulations with maximum 
distances larger than the dmax,obs then gives us the proba^ 
bility P that the observed dataset was a statistical accident. 



Because the model unimodal distributions are Monte 
Carlo realisations of the true model distributions, this pro- 
cedure (i.e. generating the 100 datasets and measuring the 
significance) was repeated for each of the 10 model grid re- 
alisations described above, to allow more accurate measure- 
ment of the significance. The significance is defined as 1 — P, 
with P defined as above. Significances for the total sample 
of 60 galaxies, and the isolated sample of 36 galaxies are 
given in Table 1. Quoted significances are the mean (along 
with the error in those means) of the 10 runs. Note that 
we also carried out this analysis in two different program- 
ming languages on two different computers using different 
sets of random number generators: the results were found to 
be indistinguishable to within the quoted errors. 
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Table 1. The significance of the bimodahty in the Ursa Major 
Cluster 



Sample 




Significance 


total 


60 


58% ± 4% 


isolated 


36 


96.0% ± 0.5% 


isolated (re-fit) 


36 


86% ± 1% 


low i isolated 


23 


42% ± 3% 
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Figure 2. Uncertainties in inclination corrected surface bright- 
ness against the axial ratio b/a. We show the surface brightness 
of the galaxy corrected for inclination assuming an infinitely thin 
disc (solid circles) connected by a dashed line to the inclination 
corrected surface brightness assuming an intrinsic disc axial ra- 
tio go of 0.20 (open circles). Error bars denote the uncertainty 
in surface brightness due to (a rather generous) error in axial ra- 
tio of 0.05. Dotted lines connect the inclination corrected surface 
brightness to the uncorrected surface brightness. 

4 DISCUSSION 

Prom Table 1, it is clear that the bimodality shown by the 
total sample is not statistically significant. However, the iso- 
lated subsamplo shows a bimodality significant at the 96 per 
cent level. To understand if this high significance of the his- 
togram for the isolated subsample reflects a true, astrophys- 
ical bimodality in the Ursa Major Cluster, it is necessary to 
understand the assumptions used to construct the isolated 
galaxy SBD. 

An obvious check is to sec whether systematic effects 
in the data can produce a bimodality, e.g., because of bulge 
contamination in the fits. We re-fit the surface brightness 

§ Note that the binning of the data in a histogram also, strictly 
speaking, invalidates the use of a standard KS test, although if 
the binning is over a much smaller scale than the 'structure' in 
the SBD (as it is in this case), tests have shown that the effects 
on the final significances are negligible. 



profiles presented in TuUy et al. (1996) using the 'marking 
the disc' technique (the same technique that was used by 
TV) , to check that there were no obvious personal biases in 
the fitting of the surface brightness profiles. Comparing our 
K' band central surface brightnesses with those of TuUy et 
al. (1996), we found a mean offset of +0.16 magarcsec"^, 
and that 68 per cent of the central surface brightnesses com- 
pared to better than ± 0.25 magarcsec"^ (both with and 
without this mean offset applied). When corrected for in- 
clination in exactly the same way as TV (assuming an in- 
finitely thin disc; by subtracting 2.5 logio(6/a), where b/a 
is the minor to major axis ratio), the two versions of the 
isolated subsample SBD arc virtually indistinguishable (sec 
Fig. 1, panel c), and the significance of the re- fitted distri- 
bution is hardly different, decreasing slightly to 86 per cent. 
This suggests that there are no systematic personal biases 
in the fits TV use that might introduce a spurious bimodal- 
ity. Note that dc Jong (1996a) shows that the 'marking the 
disc' technique is unbiased, with --^ 90 per cent of the central 
surface brightnesses accurate to 0.5 mag, compared to more 
sophisticated bulge-disc decomposition techniques. 

However, because there was no inclination cut applied 
to the Ursa Major sample, some of the inclination correc- 
tions that were applied to the K' band surface brightnesses 
were rather large: in ~ 1/3 of the cases (in either the total 
or isolated samples) they exceeded 1 mag. However, large 
inclination corrections are highly uncertain for a number of 
reasons, even in K' band. 

In Fig. 2, we show only two of the many uncertainties 
that affect inclination corrections in K' band. First we com- 
puted the corrected surface brightnesses, assuming i) an in- 
finitely thin disc (solid circles) and ii) an intrinsic disc axial 
ratio qo = 0.2 (Holmberg 1958; open circles). These two val- 
ues are connected with a dashed line in Fig. 2. The error bars 
in Fig. 2 represent the random uncertainty in the inclination 
correction due to axial ratio errors of 0.05. Taken together, 
this shows that the inclination corrections for galaxies with 
axial ratios b/a below 0.4 have considerable systematic un- 
certainties stemming from our poor knowledge of the intrin- 
sic axial ratios of galaxies, coupled with the random uncer- 
tainty associated with measurement error. Taken together, 
consideration of these two uncertainties argues that the in- 
clination corrections applied in this paper and by TV are 
likely to be lower limits. 

However, a number of effects conspire to make inclina- 
tion corrections smaller, suggesting that the inclination cor- 
rections might actually be upper limits. Firstly, the effects 
of dust in nearly edge-on galaxies is non-negligible, even in 
the K' band, which could lower the inclination correction 
by roughly as much as 0.3 mag (TuUy et al. 1998). Secondly, 
averaging of the surface brightness profiles over elliptical an- 
nuli at high inclinations (as in this case) is likely to produce 
systematic underestimates of the real central surface bright- 
ness (Huizinga 1994). Finally, at large inclinations, the as- 
sumption of a thin, uniform slab disc breaks down: vertical 
and radial structure in the galaxy will affect the inclination 
correction in non-trivial ways. Accordingly, in Fig. 2 we have 
also connected the inclination corrected surface brightness 
with the uncorrected surface brightness with a dotted line. 

In Fig. 2, we can see that uncertainties in these large 
inclination corrections may go some way towards filling in 
the gap in the SBD. As an example, we consider the four 
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galaxies with the largest inclination corrections. UGC 6667, 
UGC 6894, NGC 4010 and NGC 4183 have > 2 mag in- 
clination corrections and uncorrected (reasonably high) K' 
surface brightnesses between 17 and 18.5 magarcsec"^. The 
largo corrections move thorn across the gap from the high 
surface brightness into the low surface brightness regime. 
However, if these corrections are at all uncertain, they may 
in reality lie in the gap of the K' band SBD. Moving even 
one or two of them into the gap in the isolated galaxy SBD 
would docroaso markedly the significance of the bimodality 
in the isolated subsample. Certainly a few of the galax;ies 
mentioned above suffer from problems that may seriously 
affect the value of the inclination correction (not because of 
the uncertainty in axial ratio, but because of uuccrtaiiities 
in e.g. the vertical and radial structure in the galaxy): NGC 
4010 is clearly lopsided, with one half of the galaxy puffed 
up with respect to the other half, while NGC 4183 has a 
clear warp. 

Because of these serious uncertainties in the larger incli- 
nation corrections, and because there are no clear recipes for 
producing more realistic inclination corrections which take 
into account dust and the vertical and radial structure of the 
disc, the only fair thing to do is omit the high-inclination 
galaxies from the SBD. Accordingly, the high- inclination 
galaxies were removed from the isolated subsample, leaving 
23 galaxies with b/a larger than 0.4, corresponding to an in- 
clination of less than 66° or inclination corrections smaller 
than 1 mag. The resulting SBD is shown in Fig. 1, panel d: 
the bimodality is insignificant in this SBD, due primarily to 
small number statistics (Table 1). To check how large the 
sample size needs to be to detect significant bimodality in 
the low inclination isolated subsample, we simulated intrin- 
sically bimodal datasots and tested them using the above 
procedure. Bimodal datasets were generated using two nor- 
mal distributions with a = 0.5 mag and means of 17.15 and 
19.85 magarcscc"^, with the high surface brightness peak 
containing on average 63 per cent of the galaxies (thus re- 
sembling the TV distribution). These simulations showed 
that in order to detect a bimodality at the 95 (99) per cent 
level in the low inclination sample at least 4/5 of the time, 
the sample size needs to be increased to around 50 (100) 
galaxies. 



5 CONCLUSIONS 

We have re-analysed TuUy & Verheijen's (1997) bimodal 
SBD, using a Kolmogorov-Smirnov style of test to estimate 
the likelihood that a single-peaked distribution would be 
able to reproduce the bimodality in the Ursa Major Cluster. 

• The total sample of 60 galaxies is inconsistent with a 
single-peaked distribution at only the 58 per cent level. 

• However, the isolated subsample of 36 galaxies is signif- 
icantly bimodal: it is inconsistent with the null distribution 
at the 96 per cent level. 

• We re-fit the K' band surface brightness profiles of the 
Ursa Major sample, and found that this re-fit made relsr 
tivcly little difi^crcnco, with the isolated subsample retaining 
a reasonable significance of 86 per cent. 

• However, we argue that large inclination corrections are 
uncertain even in the near-infrared, placing the reality of 



the gap between high and low surface brigfitness galaxies in 
some doubt. 

• When the galaxies with inclination corrections > 1 mag 
are removed from the isolated subsample, the significance of 
the 23 galaxy subsample drops to only 40 per cent. 

• Assuming that the low inclination, isolated galaxy SBD 
is truly bimodal, in order to increase the significance to 95 
(99) per cent, it is necessary to increase the low inclination 
isolated galaxy sample size by around a factor of two (four) . 

To summarise, if inclination corrections greater than 1 mag 
can be trusted, then there is reasonable evidence that the 
Ursa Major Cluster SBD is bimodal. If, however, these large 
inclination corrections are in some doubt, then the SBD 
lacks a sufficient number of galaxies to argue convincingly 
against a unimodal SBD. Either way, to convincingly demon- 
strate at the 99 per cent level that the near-infrared spiral 
galaxy SBD is bimodal will require an independent dataset 
to be used with around a factor of four more low inclinar 
tion, isolated galaxies than are included in the Ursa Major 
Cluster sample. 
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